3 Relativistic wave equations 

3.1 Klein-Gordon 

The Schrodinger equation 

= Hi; (3.1) 

2 

H=^ + V (3.2) 

contains first derivatives with respect to time, t, and second derivatives with 
respect to (xyz). We now wish to handle all four coordinates (xiX2X3a;4) on 
equal footing. 

The Klein-Gordon equation was introduced by Fock (1926ab), Gordon 
(1926), Klein (1926-1927), Kudar (1926) and Schrodinger (1926). We wish 
to satisfy 

£;2 = cV + mV (3.3) 

(Compare with (1.40), m = mo). 

Taking the square root, we would get 

H = VcVTmV + V (3.4) 

but the problem would be to actually take the square root. 
Reorganizing, we get 

{E - F)V' = VcV + (3.5) 

Squaring (operating twice) 

{E - Vfi/j = (E^ - 2EV + y2)V' = (cV + m^c^)V; | : 2mc^ (3.6) 

+ ^"^'^ - 2^^^' - 2EV + V')]i. = 0\E = mc' + e 

\^ + \m(? - — "^( m^c^ + 2tm(? + - 2mcV - 2tV + F^)] ■0 = 
'-2m 2 2mc^ ^=^= ^=^= 

^ V ' 

hm (mass- velocity) of Pauli 

hm is often expressed as 

4 

iT"m 



The external electromagnetic field A,V can be included through the 
minimal substitution p ^ p — ^A, q being the charge of the particle. For 
e = |e|, (3.6) becomes, for electrons, 

e 
c~ 



[E - Vfi) = [c^{p + -Af + rr?c^] i/j (3.8) 
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[This is valid in the Gauss-cgs system of units or in atomic units. In SI, use 

p-^p-qA] 

For a free particle, replacing 

(3.9) 

and introducing the d'Alembert operator 

° = V^-^£' (3-10) 



c 



& 



,2 



r - ^2 — - rri^c^ + h^c^V^] ^ = I : h^c^ 



[□-^]^ = (3.11) 

The equation (3.7) is good for a pionic atom, se e.g. K-C Wang et 
al. PRAS^(1980), 1072. It gives a wrong answer for the fine structure of 
electronic atoms^. 

Multiply now (3.11) from left by V'*: 

^*[U-'-^]iP =0, and I + 
=0 I - 

V'*(nV') - (□V'*)^ = (3.12) 



Using the definition of □ 



^*(^' - ;^^)^ - [(^' - ^2^2r]i> = (3.13) 



As 

rm - (r)> = [ri^' - (r)V]' (3.14) 

and 

V'*(VV) - V'(VV*) = V • - V'(Vi/'*)] (3.15) 

we get the continuity equation 

^ + V-l = (3.16) 



where 



3:{r,t) = -^(V;*VV- (V^*)V) (3.18) 



'The proof is left as an exercise for the reader 
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reduce to the right non-relativistic hmit (see e.g. Kvantkemi I, p. 20-21). 
Now the p is still real, but not necessarily positive definite. Replace in (3.17) 

E 

p= — ^|V^|2. Now, if £; < 0, p < (3.19) 



3.2 Dirac 

We are searching for an equation of the form 



(3.20) 



Because of the first derivative ^ , we would like to have first derivatives ^ 



etc. as well. Note that the series expansion of the square root in 



{E - V)i; = mcVl+P7c2 

would contain all powers of {p^ /(?)■ 

Suppose now that the Hamiltonian is linear in all and that the wave 
function * has N components, 



Then the most general free-particle wave equation is 

d 



(3.21) 



(3.22) 



where n runs from 1 . . . iV, and 

..f, . d . . d . . d 

In terms of the N x N matrices a and /3, 

d 

ih—"^ = \-ihca ■ V + /?mc^l* = Hd"^ 
at 



with the Dirac Hamiltonian 



ho = ca - p + I3mc 



(3.23) 



(3.24) 



(3.25) 



The D in hu stands for Dirac. 

The components of a are the N x N matrices (A^ > 4, see below) ax, oty 
and az- In order for ho to be Hermitian, a and (3 must be Hermitian: 



= a, 



(3^ = 13 



(3.26) 
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For all points in space-time to be equivalent, a and /3 must be constant 
and dimensionless. Consequently they commute with r and p. 
We still want to satisfy 



for all components jpi . . .i/jn- 



-h"^^ = [-n^c^V^ + m^"]^ (3.27) 



Then the Dirac equation 

ho^ = ih—'^ 
at 

will connect the different components while every individual component, ipi, 
will satisfy the Klein-Gordon (K-G) equation (3.27). 



3.2.1 Why must AT > 4? 

Starting from 



Pai = -ail3 1/3 x 
ai = —j3 aiP (as = I) 

A B 



we see that 



Tr AB = Tr BA 
TV = -TV {aiP)(3 \f = I 
= -TV ai 
= 

The Ui were ±1. We have the same number of both =^ N must be even. On 
page 17 it was shown that 2 was too small ^ N >A. □ 

3.2.2 Properties of Dirac matrices 

Multiply now 

(m- - ho)^ = 

from the left by the operator ih^ + ho ^ 
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in other words 



q2 



-i^c E ti + otiff) £: + m^cV'] * (3.28) 



where we have used the fact that ^ = 0. 
Prom (3.27), 



' aiOj + ajUi = 26ijl 

Pai + aiP = (3.29) 
/?2 = 1 

I being the N x N unit matrix. Because a and P were Hermitian, their 
eigenvalues are real. According to (3.29), the squares of these eigenvalues 
equal 1. Hence the eigenvalues are ±1. Dirac noticed that (3.29) is satisfied 

by 



P=(^ J ) (3.30) 
a 



(3-31) 



where 



'H-)' -H-)' "-(""o)' <^-^^) 

3.3 Some properties of the Dirac equation 
3.3.1 Solutions with ±E 

For a free particle, 

[hD,p\=0 (3.33) 
and hence the two operators have simultaneous eigenfunctions '^{p,x): 

ho^ = E{p)^ (3.34) 

Operating again with ho, 

hjj^ = E{pf-^ (3.35) 
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with the solutions 



(3.36) 



corresponding to " electron- Uke" and "positron-like" solutions, respectively. 
If only one of them existed, ho would be a multiple of I, and commute with 
all matrices. 
As (3.29) ^ 

{a, 13} = ^ 

[a, /?] =al3-l3a = -2/3a ^ 
[ho, (3] = -2cl3a-p^0^ 
3/3 so that [hd, (3] ^ ho 7^ aI^3+E and -E 

3.3.2 Inclusion of electromagnetic fields 

Use again the minimal substitution 



[In Gauss-cgs. In SI, use p — eA. Now e < 0.] The non-relativistic "dia- 
magnetic" term is recovered here in 2'^'^-ordcr perturbation theory with in- 
termediate positron-like states. Using the closure. 




(3.37) 




(3.38) 



n 



AE = 2mc^ 



we get 



n 



= e'^{ip+\A^\^+)/2mc^ (3.39) 
the usual diamagnetic term. (Recall that ^{p— ^A)'^ the term g^^J;^.) 



3.3.3 Free-particle solutions 



In atomic units (e = nie 



h = 47760 = 1, c = 137.036), 



ho = ca- p + (?f3 



(3.40) 



Using the Ansatz 



= u{p)e'P-^ 



(3.41) 
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the Dirac equation ho^ = gives 



-E + c^ 







CPz 



-E + c^ c{px + ipy) 

cpz c{px - ipy) -E - (? 
c{Px + Wy) -cPz 



- iPy)' 

-CPz 



-E-c2 





Ui 












U4 



(3.42) 



This hnear, homogenous equation has solutions if det [] = 

{E''-c^-cVf = 

E = ±c^{p^ + c2) = ±Ep 
Denoting Ep + by A, the spinors become {az= the spin): 



(3.43) 
(3.44) 



{E,az 



{-Ep, i) 



( -^P' 2' 



Ui 



ui 1 

U2 1 

tt3 cp^/^ c{p.^ - ipy)/A 

U4: c{px + Wy)/A -CPz/A 



-CPz/A -C{px - Wy)/A 

-c{px + iPy)M CPz/A 



1 







1 



The normahsed spinors become 

u{p,Ep) -- 
u{p, -Ep) = 



2 En 



1-Q ^ 



2Ep VV': 



/-^M^(a\ 



^2 



(3.45) 

(3.46) 
(3.47) 



where (^) is cither (q) or Cj*). 

The ratio of the norms for £'p > 0, 



11)1 {E + c^f (2c2)2 4^mc^ 4^c 



(3.48) 



for small v/c. Therefore ipi and ip2 are called the large and small compo- 
nents, respectively. 

For electrons in light atoms, v 1 a.u. and this plane-wave estimate 
gives 



A^c' 4^37 
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3.3.4 Probability density 

For satisfying 

d 

ih— = -ihca ■ V* + PmcH 1*^ x (3.49) 
at 

where the Hermitian conjugate 

*t = (V;*,^*,V'3,V'a) (3.50) 

satisfies 

d 

-ih—^^ = ihc{V^^) ■ a + mcH^fi \ x * (3.51) 
Performing the two multiplications and subtracting, we get 

ih— (*t*) = -inc(*^a • V* + V*^ • a*) (3.52) 
1 d 

- — (*%) + V- (^^a^') = (3.53) 
cot 

Comparing with the classical continuity equation 

^ + V-i = (3.54) 



we identify 

p = q>^^ = (^t^i + ^t^2 + • • •) (3.55) 
I = c^^a^ (3.56) 

Prom j = pv, the velocity operator (local velocity of the Dirac electron) 

v = ca (3.57) 

3.4 The Pauli limit 

By substituting E = mc? + £, the Dirac equation becomes 

[ca-p+ {(3 -I) mc^ + V]^ = (3.58) 


^0 -21 

In terms of the two-spinors 



cg_ ■ pip2 + V'ipi = ^V'l 
ca ■ pipi + (— 2mc^ + V)tp2 = £^'2 
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(3.59) 



-02 



ca ■ p 

2mc^ + £ — V 
q_ ■ p 



2mc 



(3.60) 
(3.61) 



(Note that (3.60) is still exact.) Using the above approximation with tpi 

tpnr, (V'llV'l) = 1 ^ 



{^\^) = (V'l|^l) + (i/'2|V'2) 

= 1 + (V'2|V'2) \p-a^a-p = p^l 
1 



1 + 



(V'lblVi) 



4m 2 

Then the expectation value oi ho becomes 

{tjji\ca ■ plV'2) + ('i/'2|co: • pIV'i) — 2mc^('i/'2|V'2) 

+ (V'll^lV'l) + (^2|F|V'2)]/[l + (^21^2)" 
r 1 1 Imr^ 



(3.62) 



+ {tPi\V\'il^i)+ ^ 



r(^i|a-pFa-p|V'i) /(*|*) 



1 - 



1 



1 -1 
{ipi\g_ ■ pVa ■ pli^i) 

(?/'i|p^|V'i) + 0(c""^) l/c"^ and higher order terms neglected 



1 



4m^c^ 
Trick 1: 



(3.63) 



p^ 



F^l = (.nr - ^)Vl 



(3.64) 



S — Snr 



r£nr(V'lb^lV'l) 



4m^c^ 



(3.65) 
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Trick 2 : 

GiGj = i(5yfe(7fe + Sij, (3.66) 
a ■ pVa ■ p = aiPiVajPj 

= ^SijkakPiVpj + pkVpk 

= 'Ajk(^k [{Piy)P3 + VPiPj ] +PkVpk 

=0 

= fuj-(VV xp)+pkVpk (3.67) 

[Note: piPj = as i and j run over all combinations (summation); pipj — 

PjPi = 0.] 
Here 

PkVpk = ^ [pkVpk + PkVpk] 

= I [V +P^V] + ^ [{pkV)pk - PkiPkV)] (3.68) 
where the first term becomes 

liHVp^ = £nr(V'l|pVl> - ^(^ll^'^IV'l) (3-69) 

and the second term becomes {dk = ^) 
^ (^1 1 (Pk V)pk - PkiPkV)\tpi) 



^{i^i\{pkV)\pkipi) - ^{pk^i\{PkV)\ipi) 

' 2 
" 2 



Pk = -ihdk 



[ J ri{dkV){dki^i) + J (5fe^i)*(5feF)Vi] 

{[riiOkV)^!]":^^- Jidkri){dkV)i^i- J rMv)^! 



+ 

1 



J (9feV'l)*(S/c^)Vl} 



= ^^(00)^0 (3.70) 

Assembling the pieces in (3.65), (3.67), (3.69) and (3.70): 

1 r, 

+ - 



4^2c2L^V'i|^^-(Vl^X|,)|V^i) 

£nr + (/im + hso + /id) + (^(c"^) (3.71) 
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' {hm) = - ^\,^2 (V'll/lV^i) (mass-velocity) 

^^so) = ^^2g2 ^^i|g ■ ^ (spin-orbit) 

(^d>= ^^^(V'llWlVi) (Darwin) 
Returning to (3.65), in terms of il)2, 

£ = £nr-£nr + (V'2|1^|^2) 

For the Coulomb potential of a point charge Z, 

VV = -ZV2(^) = 47rZ ()(r) ^ 

3-D <5-func 

(^d) = ^(V'i|(47rZ)5(r)|V'i) 



Example 1: For the Is state of a hydrogenic atom, in a.u.. 



TT 



- ^ {{E-Vf) 



1 



2mc^ 



(f )^-2(-f )(F) + (F2) 



As 



^-4Z^^^-Vd.^(-4Z^)^^-Z^ 
Jo \r y 2Z 



(ftso) = 



' 1 + 2 



8c2' 



31 



Summing up, 



5 4 



(3.76) 



3.5 Central fields 

We now return to the full Dirac equation with a local potential V{r) = V{r). 
As 



[l, hi)] = ica X p 
[s, /id] = — ica X p ^ 



(3.77) 
(3.78) 



neither Z or s can have simultaneous eigenstates with /id but their sum 

1 = 1 + s (3.79) 

can because 



[l,hB\=0 



(3.80) 



Details: 



I = rx p 

[/,/?] =0, [l,V]=0, [l,a-p\^0 



Consider the component h: 



[li,a-p]= ai [h,pi] +a2[li,P2] + asihyPs] 

=0 

i J k 



l = rxp = ~i 



X y z 
dx dy 



h = -Kydz - zdy) 



[h,P2\ = {-if[ydz - zdy,dy] 



92 52 „ 92 q2 



( ^ ^ ^ 

\ ovoz ov^ 



dydz 5y2 dydz dy^ 

d . ~ 

— = ip3, cycl. 
oz 



[h,QL-p\ = a2ip3 - "31^2 = i(a X p)i, cycl. 

\l,a - p] = \ ax p □ 
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Similarly for [a, a ■ p] (=with spin), 
[ai,a-p] -- 



As [(71, (T2] = 2i(73, cycl.. 



(71 \ / o:-p\ _ I (7 • p\ /5i 
aij \a-p J \a-p 7 V 

aia ■ p\ fO a - pai 
aia ■ p J \ pen 

[cri,o:-p] 
[(71, a -p] 



[(71, a • = (7l((7ipi + (72P2 + (73^3) 

-((7lPl + (72P2 + (73P3)(7i 
= Pl + (7l(72P2 + (7i(73P3 
-Pl - (72(7lP2 - (73(71^3 

= |_'7i^^P2 - ba^^xl P3 

2i<T3 2icr2 

= 2i[(73P2 - cr2P3] = 2i(p2(73 - P3C72) 
= 2i(p X a), cycl. 
[ct, a- p\ = 2ip X a = — 2ia x p □ 



Define 



(3.79) 



s = ]^ha (3.81) 



As a further detail, consider the operator 

K = {Va ■1 + 1)13 (3.82) 

We shall show that 

[K,j] = (3.83) 

[K, hr>] = (3.84) 

In other words, K, /id, and jz can have simultaneous eigenstates. 
Theorem a: 

For two vectors A and B which commute with a but not necessarily with 
each other, 

{a- A){a- B) = A- B + ia- {Ax B) (3.85) 
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Proof : 

The a matrices simultaneously satisfy 

af = 1 {i = k) 

a^ak = iSikm {i / k) 



[o:,A]=o 



= GiGi AiBi + ^ aiajAiBj 



iSijkSk (4xB)fc 

= A - B + ia - Ax B □ 
The equation holds for both two- and four-component a 

Theorem 6: 

{a-A){a-B) = ia-A){a-B) 

Proof : 

/ a-A\( a- B\ _ (a- Aa- B 
\a-A ) \o:-B /V o:-Aa-B 

= {a^-A){a^.B) □ 

denotes a four dimensional a. 
Theorem c: 

(a • • = -75^ • B ia • (A X S), 

where 

75 = 71727374- 
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Proof: 

We use the representation 



Then, recalhng that i' 



,•3 



Then 



75 







7 



75 



-CTl 





P 





0-2 



-ia 
ia 

1 
-1 



-(72 





-cri(T2 






(71(72 

(71(72(73 
(71(72(73 

1 

1 



75" 



"75 



75a = - 



4 

q: 75 



1 

1 

a 
a 

1 

1 

a 

(7 



a 
a 

1 

1 

a 
a 

1 

1 





(73 



-(73 




1 
-1 



(73 
(73 



(71(72 = 1(73 



a 
a 

a 
a 

a 
a 

a 
a 



= —(7 



= —(7 



= —a 



—a. 



Multiplying Theorem a, (3.85), from left by —75: 



75(7^ ■£) {o-B) = -75 A • B + i(-75(7^) • A X ^ □ 



Theorem d: 

Because 

I = rx p (-L p), 
p.l = l.p = 0. 

Furthermore, 

lxp + pxl = 2ip 
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Proof : 



Ax B 



Ax B + B X A 



i J k 

Ax Ay Az 

Bx By Bz 

l{[Ay,Bz\-[Az,By\) 

+j{-[Az,Bx] + [Ax,Bz]) 

+ k{\Ax,By]-[Ay,Bx\) 



{Ixp + pX l)x = [ly:Pz] " [lz,Py] 



d_ 

dx 
d 



d d 

Then the anticommutator 

{a-l_,a-p} = {g:- l_){a- p) + {a- p){a- 1) 

(3-90) r, „ , „ 7 . : 4 



75[1 • P + P • 1 ■ {Ix P + px I)] 

=0 =0 2ip 

2a ■ p (3.97) 



Theorem e: 



Because 



[K, a- p\ = [Pa ■ l_,a- p] + [P,a- p] 



—a[3 = (3a, 

-a ■ pPa ■ I = P{a ■ p)a ■ I =^ 

[K,a-p\ = P{a-l,a-p} +2(3a • p = 
V ' 

—2a-p 



(3.98) 



Theorem /: 

As 



[/3,/?] = 0, [i^,/?]=0. 



(3.99) 
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Theorem e+/ 

[K,hB] = (3.100) 
It remains to show that [K,j] = 0. 

Theorem g: 

Prom theorem a, (3.85), taking B as an arbitrary vector, commuting 
with A, 

{oi- A)a = A + iax A, (3.101) 

where we used 

a - Ax B = ax A - B. 
Similarly, letting A be arbitrary, 

oiia- B) = B + iBx a. (3.102) 

Letting then A = B = IsLud subtracting (3.101)-(3.102), 

■ L ^] = i(TX I — ilx a 

= 2iaxl (3.103) 

Theorem h: 

Let A commute with I. Then 

3 3 

\ A • /, Z] = ^ ^ ''^^'[AjliljXj ljXjAjlj\ 
i=i j=i 

ij "^"^ — — ' 

= -i4 X L 

Letting A = a, 

-1,1] = — io: X / 

Theorem i: 

The desired commutator (note: [K, /3] = 0) 
[K,i] = [Pa-l,l + ^a] 

= K-io: xl) + ^P{2ia X /) = □ (3.106) 

Summarizing, 



(3.104) 
(3.105) 
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[K,j] = [K,hB] = \j,hB] = 



(3.107) 



Denote the eigenvalues oi K by — k: 

Kx": = -/^X^, (3.108) 
where the a (in K) are 2-component ones, and the 2-component spinors 



X« = \lsjm) 
The quantum number k carries both j and I: 



K = < 

The lowest values are 
K = 



-l-l, j = 1+^ (state T) 

/, j = l-^ (state T' or 'r 



Prom 



-1 1 -2 2 -3 3 -4 

Sl/2 Pl/2 P3/2 C?3/2 C?5/2 /5/2 /7/2 

s p* p d* d f* f 



j = L + s, 



•2 



1^ + 21 • s + 



2 s^] 



Consequently 



K = a-l + l 

= K\lsjm) = [jU + 1) - l{l + 1) - + 1) +l]x™ 



whence 



(3.109) 



(3.110) 



(3.111) 



(3.112) 



K = Z(Z + l)-(i + -)2, 



(3.113) 



as can be verified from Table (3.111). 
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For the 4-component case we introduce the Ansatz 

* = = ("^A = ( '^'^^^ ) = '-( ) (3 114) 



for which, indeed 



^^:r = (-'n^' ^(t!)=-^^'^■ (3-115) 



-a • i - ly VV'2 
This * is an eigenfunction to hD,K,j'^ and jz- 

3.5.1 The radial Part 

From 

Ax{BxC) = {A- C)B - {A ■ B)C, 

{A- B)C = {A-C)B- Ax {Bx C), (3.116) 

with 

A = B = r = {i,j,k), C = V, 
V = f (r • V) — f X (f X V) 

= er^-i-xl. (3.117) 
or r 

Here Cr is the unit vector along the radius (and I = rx p). Consequently 
the kinetic energy term 

d 1 

a-p = —\arT, a-rxl. (3.118) 

— or r 



Denoting 

V 
,0 a 



and using 

a- Aa- B = A- B + \a- Ax B [3.88-3.89] 

for A = f and B = I, 

ia-r){a-l) = r^+ia-rxl (3.119) 

ar =0 

(f • [ = as the angular momentum has no radial components.) Multiplying 
by i75 and recalling that a = — 750:, 

i75CKr(« • I) = —75^ • rx I 

i75(7r(a-i) = a-rxl_. (3.120) 
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Thus the Dirac Hamiltonian 



h-D = ca - p + f^m.c^ + V{r) 



iC75crr 



d__l 
dr r 



a-l 



+ I3mc^ + V{r). 



with 



K = ^{a-l + l), 0^ = 1 
a-l = l3K-l, 



hi) = icj^ar 



9,1 f^K 
dr r r 



+ Pmc^ + V{r) 



Theorem a: 



al = 1 



(Uke in the cartesian case, af = 1) 
Proof : 



(Jr = a; • f = ^ ai^{axX + ayy + a^z) 

i 

al = fraclr^ + ^ XiXj {ajaj 



i>j 



= 1 



Theorem 6: 



{K,ar} = 



Proof: 



K = f3{a-l + l) 
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(3.85) 



{ql ■ r){a ■ 1) = r ■ I + ia ■ {r X I) 




{K,ar} = {P{a-l + l),a-f} 

= Pia- {Ix r + rx I) +2/3a • r 



= □ 



Theorem c: 

CTrX"^ = -X"^. (3.127) 



Proof: 



= l^{(^rXK) 
Thus the eigenvalue of (TrX^ is while 

KXn = -/^XlT 

whence 



C7rX'^ = -X-. □ 



Recalling now that 



75 



1 

1 



(3.124) becomes 



ic(-) 



1 

1 



d_ 1_1 A 
(?r r r VO — 1 



1 
-1 



1 
-1 



(3.128) 



[3.108] 



mc^+V-E 



9X^ 



0. 



where K'^ denotes the four-dimensional K-operator. Recalling 



4/ 9Xk 



(-K) 



f 9Xk 

VfX-K 
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d_ V 



+mc 



.9X-K 
9Xk 



9Xk 



+ {V-E)[^ 



For the two components 



+ 



i/X-K 



{±}mc^ + V -E 



0. 



i/X-K 



Dividing by the angular part and i we get 



-c(— + 
\dr r 



+ 



1 + K 



Introducing 



mc^ + V - E)g =0 
g + {-mc^ + V -E)f =0 

P = 

Q = r/' 



dP _ dg 
dr dr 



dQ 
dr 



f + r 



dr 



f dy _ 1 dP _ P 

dr r dr r^ 
d/ ^ IdQ _ Q 

dr r dr r^ 



--j I-(I-k)— + -mc H — = 

rVdrr r J \ c J r 

1/dP P , / ^ 

- -5 + 1 + K - + (-mc+ ) - = 

^rVdr r r J \ c J r 
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dP P f V-E\^ 

-;- + K- + -mc+ ]Q = 

dr r \ c 



— K h — mc H P = 

^ dr r \ c J 



E-V^ 



3.5.2 Non-relativistic limit 

Let the potential \ V\ <^ E ^ mc^. For the E > states 



-mc + 





— 2mc =^ 


c 






1 /dP 




2mc V dr 


dQ _ 


1 (p. 


dr 


2mc V 



r / 



J 

r 



Denoting = mc^ + £, for bound states e < 0, 

_r r^ r \ r / 

mc^ + £ — F 



+ (2mc) ^— mc + 



P" - ^ ' P + 2m(e - y)P = 0. 

-^-1 . X fH-i)H) 



P + 2m(£ - V)P = 0, 



the usual non-relativistic equation. 

3.6 The Dirac-Coulomb Problem 

(Darwin 1928, Gordon 1928). Consider 

Z 



V 



For normalized bound states 



/ r'^dr{f + g'^) = / dr(p2 + Q^) = i 
Jo Jo 
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(3.134) 



dP P f Zl 

— + K— + -mc Q = 

dr r \ c r c / 



dQ Q 
dr 



r V 



-mc + 



^1 ^ 
+ - P = 



(3.142) 



c r 



So far /i = 1. Let now m = c = 1, as well. Then the unit charge 
e = a = 1/137.036. Furthermore denote 



\ = ^{l- E") 
X = 2Xr 



(3.143) 
(3.144) 



and introduce the new functions (pi and 02 for which 



P = VlTEe-^'-{(l)i + ^ 



(3.145) 



Then 




-A(<^1 + (f>2) + 
-A((/.i - (1)2) + 



#1 


dr - 


dr 


d0i 


d(p2' 


dr 


dr - 



(3.146) 



Straightforward algebra on Eq. (3.142) gives: 



Vr+^e-^^[-A(<^, + <^2) + ^ + ^ 
L dr dr 



_^!VrqrEe-^'-(0i + (/.2) + 

r 



£; + - + i 

r 



Vr^e-^'-(0i-(/>2) 



Dividing the top equation by Vl + Ee and the bottom one by ^/l — Ee 
and noting that ^ = 2A^ we get 



-A(0, + <^2) + 2A(^;fi + ^;^ 
dx dx 



('^1+02) + 



2\Z 

E+ + 1 



X 



H - 02) - 



■ 2XZ 
E+ 1 



X 



1 - 

l + E 



l + E 
1 - E 



(01 - 02) 



(01 + 
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Dividing now by 2A, gives 



1 K 

4>1 + ^'2 =7){<l>-i-+ 4>2) (01 + 4>2) 



+ 



■ 2XZ 
E+ + 1 



2^{1 + E){1- E) V 1 + E 



1-E 



1 K 

^'1 - (p2 =2^^^ ~ '^2) + -i4>i - 4>2) 
E^ 1 



2^{l + E){l-E) \ l-E 



l + E 



{(pi - 02) 



(01 + 02) 



Now, remembering the definition of A, (3.143), the above can be reorganized 
to 



01 + 02 



1 K E+^+1 

2^ 2{1 + E) 



+ 02 



1 K E+^^ + 1 

2~ x' 



2{l + E) 

2\Z 2 



2 X 2{1-E) 

+ 02 



1 K 



2 X 



2(1 - E) 



Writing down the sum and difference of the above { (sum at top, difference 
at bottom): 



20'^ =01 

+ 02 
202 =01 

+ 02 



1+ {E + 



2\Z\ 1 / 1 



X 



nil 



1 \ 1/ 1 1 



2k 



X 



+ 



(^E + 



+ E 1-EJ 2\1 + E 1-E 

1 / 1 



2AZ\-1/ 1 1 
+ 



2k 

X 



X J 2 \l + E 1-EJ 2\1 + E 

/ 2AZn1/ 1 1 N 1/ 1 _ 1 \ 

"^^J2Vl + £;'^ i-e)'^2\i + e~i-e) 

V X J2\l + E 1-EJ 2\1 + E 1- 



-1 

1 + E 1-EJ 



As 



-2E -2E 



1 + E 1-E l-£;2 A2 
1 1 _ 2 _ 2 

1 + E ^ 1-E ~ 1-E^ ~ A2' 



E. 
(3.147) 



(3.148) 
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2\Z\/-E\ 1 



+ 4>2 



2k / -i\Zj\ 

E + — 

X \ X J 



2\Z\ 1 E 
A2 + A2 



2<^2 



2n ( ^ 2\Z\ 1 E 



+ 02 



i-(e + 

\^ ZE- 
xX . 

K Z 

I X x\ 



2\Z 






X 


) 




+ 02 








,T 






\ZE^ 


" +0 


2 


- x\ - 







A2 



Introducing the series 



TO=0 

oo 



TO=0 

oo 



(3.149) 



(3.150) 



(3.151) 



m=0 



J](7 + m)«— 1 

m=0 



m=0 
oo 



xA / 



m=0 

oo 



m=0 



X .tA 



+ Pm^ 



X xA 



xA 



(3.152) 



The coefficients for every power, x'^~^"^ ^, must vanish for the expansion to 
be true: 



/ N ZE ^ / Z\ 
am[l + m) = am-i - am— PmyK + —j 



A 



^ Pmil + m) =am (^-K + j) + Pr. 



ZE 



(3.153) 



If am-1 is known, these equations gives {am,Pm) For m = (a_i = is 
assumed) , 

( / ZE\ ^ / Z\ 

/ Z\ ZE 
^ 7/?o = ao[-K+ —j + Po^- 



(3.154) 
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This linear, homogenous pair has solutions if 



7 + ZE/\ K + Z/\ ^ 
K-ZjX -f-ZE/\ 

^-[-x) -'^ +(a) =° 

72+(f)V-^^)-'^^ = 



7' = - 



(3.155) 



(3.156) 



It turn out that 



7 = +\/ — Z"^ 



(3.157) 



is the acceptable solution and 

7 = —\Jk^ — Z'^ 

the irregular one (which can't be normalized). Recall that Z included the 
proton charge Cp = a = 1/137.036. For integer Z', 



7 = +y«;2 _ (^^/)2 



(3.158) 



The Eq. (3.153b) 



—n + Z/\ K — Z/\ 

ctm m + j — ZE/X n' — m ^ 



(3.159) 



In terms of a new quantity 



, ZE 
n= — -7. 



(3.160) 



A substitution to (3.153a) gives 



ZE ( Z\K- I 



7 



Z\K-ZI\. 
m 



(3.161) 
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Here 



ZE ( Z\K-z 



Z\n-Z/\ 
m 



ZE ^2 _ 

A n' — m 

( ZE\( ZE\ o , 

(^7 + + - "I + — j + At^ - (Z/A)^ 

- (72 + 27m + m^) + {ZEjXf + ^2 _ (z/A)2 





1 




n' 




m 




1 




n' 




m 




1 




n' 




m 



-K^ + Z^ - 27m -m^ + K^ + Z^ 



E^-1 



m- 



27 + m 



n' — m 



an 



n — m 
m(27 + m 



-am-i- 



(3.162) 



(3.159) 



(3.154b) 



_ (n^-m)...(n^-l) 
""^"^ ^ m!(27 + m)...(27 + l)"° 
(l-nO(2-nO...(m-nO 
= m!(27 + l)...(27 + m) "'^^ ^^"'^^^ 

Pm ^ K- Z/X 

dm n' — m 

Pm-l K- Z/X 



am-1 n' — m + l 

I3m ??/-/» + 1 <\„i 



Pm-l n' -m am-i 

R -( \rn (n^~m+l)...(nO 
m!(27 + 1) . . . (27 + m) 



OLQ 

ao 7 — ZE/ X n' 



Z\ ^ / ZE\ 
A~J 



(-.+ f)=/3o(7 ' 



/?o -K + Z/A n-Z/X' 



(3.165) 



48 



Comparing the coefficients (3.163)-(3.164) with those of a confluent hyper- 
geometrical function, 



„ , a, a(a + 1) 

iFi{a;c;x) = 1 + -x + , -,^ 7^ H , 

c c(c + 1) 2! 

= aoa;^i^i(l - n'; 27 + 1; x) 
(/)2 = /3oa;^iFi(-n';27 + l;x) 

= ^^^ aox^iFi{-n'; 27 + 1; a?). 
n' 

For a; — 00, i-Pi(a; c; x) — e^ because 

oim -n' + m 1 

hm = nm — — >■ — . 

am-i m-»oo m(27 + m) m 

Therefore the series must terminate, That happens for (p2 with 

n' = 0,1,2,... 



(3.166) 



(3.167) 



(3.168) 



(3.169) 



(n' = is the ground state, n' = 1 the 1^* excited state, etc.) Thus n' is the 
number of nodes in ^2 (or in g) 

For n'l = 0, the coefficients of (f)i have 



(1 — n') ... (m — n') 
m!(27 + l)...(27 + m) 

1 

m!(27 + 1) . . . (27 + 



which converges. 

Introduce the principal quantum number 

n = n' + \k\. 

Then the definition of n', (3.160), gives 



"0, 



(3.170) 



ZEn 

A 



= n' + 7 



1 + 



(n' + 7)2i (n' + 7)2 



Z2 



Z2 




(3.171) 
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The K occurs only in 7 = ^/ — Z'^, +|«;| and — are degenerate for 
the same n! For instance for n = 2, 

^2si/2 = -E'2pi/2 < -E2P3/2- (3.172) 

(Note though that with QED accounted for E2s^/2 expansion 
in powers of gives (we stih have c = 1) 



_ _1Z^ Z^rS 1_ 

-t^riK "J- — ^riK — n n H 5" 7^ TTj T 

z n"^ ISn 2\k\ 
Qualitative conclusions: 



+ 0{Z^) (3.173) 



a. Largest stabilization for \k\ = 1. 

b. States with same I (and n) but different j are spin-orbit split. 

c. For |k| = 1, both / and g diverge at the origin like r'^~'^. They remain 
normalized. The isomer or isotope shifts, proportional to p{0), grow 
by a factor of 13 for U. 

d. The radial density, + Q^, has no nodes. 

e. The radial electron density suffers a relativistic contraction. 

f. Normalization is no problem 

g. Continuum, also positron-like can be solved. 

h. Integrals, properties, solvable. 

i. Literature: See Table 2.3 of RTAM I-III (1986, 1993, 2000) 

3.7 Virial theorems 

Literature: See RTAM I-III, Table 2.5. 

Fock (1930b), Gupta (1932), Rose k Welton (1952), March (1953), Kim 
(1967), McKinley (1971), M. Brack (1983). 

3.7.1 Non-relativistic case 
Let 

r' 

r^r' = Xr |r = y (3.174) 
A 

H = T + V = -^V^ + Cr" 

^ A2(^-^V'2) +CA-"(r')" (3.175) 
E{X) = A2(r)+A-"(y) (3.176) 
2A(r) -nA-"-^(F)l =0 (3.177) 



dE 



A=l 
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2{T) = n{V) 



Coulomb: n = 1 



Because E{X) has to be a minimum at A = 1, 



A=l 



> 



n{V) + n{n + l){V) >0, {V) <0 
(2 + n)n{V) > 



-2 < n < 



3.7.2 Dirac 



hi) = ca- pOPmc^ + V{r), 
V{r) = Cr"^ 



Bound states considered: 



, (0|</.) = 1, {A\A) + {B\B) 



Again, let r^r' = Ar, dr — ^ A ^dr' 



$(r) ^ A=^/2$(r') 
h-D cXa ■ p + Pmc^ + A""y(r). 

' {^'\ca-p\^') = X{^\ca-p\^) 
< ($'|/?mc2|$') = mc^{NA - Nb) 
(^'Ir-'^l^') = A-'^($|r-"|$) 



S(A) = A(ca • p) + A-"(l/) + mc^{NA - 



(-) = 







{T)=n(V) 



For n = — 1, 



(T) + (y) = c{a ■p) + {V) = 
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For n = —1, 



(y) < ^ ri(n + 1) < ^ 
-1< n < 0. 

E= {T) + {V) + mc^iNA - Nb) 
E = mc^{NA- Nb). 

E = m(?{P) 



(3.191) 
(3.192) 

(3.193) 
(3.194) 



sc- 




on D-2S 




Figure 3.1: Radial electron densities for the Is, 2s, 2p and 3s states in 
a hydrogen-like atom with Z = 80. The dotted curve denotes the non- 
relativistic density and the solid curve the relativistic density. Reproduced 
with permission from V.M. Burke and LP. Grant, Proc. Phys. Soc. (London) 
90 (1967) 297-314. 

From Figure (3.1) it is seen that the relativistic contraction is of the 
order for the states Is, 2s, 2pi/2 and 3p, whereas the contraction of 2p3/2 is 
considerably lesser. The figures 2 and 4 (inside (3.1)) show that 2s and 3s 
are procentually contracted as much as Is is. The same conclusion can be 
drawn for all s-electrons up to the valence shell in a many-electron atom. 
The p-shells also contract in the same way, but not as much. 
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